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Abstract 


Due to the exponential high gravitational red shift near the event horizon of a black hole, 
it might appears that the Hawking radiation would be highly sensitive to some unknown high 
energy physics. To study possible deviations from the Hawking’s prediction, the dispersive field 
theory models have been proposed, following the Unruh’s hydrodynamic analogue of a black hole 
radiation. In the dispersive field theory models, the dispersion relations of matter fields are modified 
at high energies, which leads to modifications of equations of motion. In this paper, we use the 
Hamilton-Jacobi method to investigate the dispersive field theory models. The preferred frame 
is the static frame of the black hole. The dispersion relation adopted agrees with the relativistic 
one at low energies but is modified near the Planck mass m p . We calculate the corrections to 
the Hawking temperature for massive and charged particles to O (nif 2 ) and massless and neutral 
particles to all orders. Our results suggest that the thermal spectrum of radiations near horizon 
is robust, e.g. corrections to the Hawking temperature are suppressed by m p . After the spectrum 
of radiations near the horizon is obtained, we use the brick wall model to compute the thermal 
entropy of a massless scalar field near the horizon of a 4D spherically symmetric black hole. We 
find that the leading term of the entropy depends on how the dispersion relations of matter fields 
are modified, while the subleading logarithmic term does not. Finally, the luminosities of black 
holes are computed by using the geometric optics approximation. 
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I. INTRODUCTION 


The classical theory of black holes predicts that anything, including light, couldn’t escape 
from the black holes. However, Stephen Hawking demonstrated that quantum effects could 
allow black holes to radiate a thermal flux of quantum particles l| . The assumption that 
the ultra-high energy modes are in their ground state was used to derive the Hawking 
radiation in the framework of quantum field theory in curved spacetime. After this discovery, 
it was realized that there was the trans-Planckian problem with the calculation ^. Due 
to the exponential high gravitational red shift near the horizon, the outgoing particles of 
the Hawking radiation originate from the extremely high (e.g., trans-Planckian) frequency 
modes. So the Hawking radiation relies on the validity of quantum field theory in curved 
spacetime to arbitrary high energies. On the other hand, quantum field theory is considered 
more like an effective field theory of an underlying theory whose nature remains unknown Q] . 
This observation poses the question of whether any unknown physics at the Planck scale 
could strongly influence the Hawking radiation. It is believed that the trans-Planckian 
physics manifests itself in certain modifications of the existing models. Thus, even though 
a complete theory of quantum gravity is not yet available, we can use a “bottom-to-top 
approach” to probe the possible effects of quantum gravity on our current theories and 


experiments 


4], One possible way of how such an approach works is via modifications of 


energy-momentum dispersion relation 

2 t-i2 2 

p = E — m , 

which has been reviewed in the framework of Lorentz violating theories in 
To study the trans-Planckian problem, the dispersive field theory models 


( 1 ) 

, 6 ] 

17] have been 


proposed, which focused on studying the effect on the Hawking radiation due to modifications 
of the dispersion relations of matter fields at high energies. These models were motivated 
by a hydrodynamic analogue of a black hole radiation 7]. Similar to the original method 
for deriving the Hawking radiation, the energy fluxes for outgoing radiation were usually 
obtained by calculating the Bogoliubov transformations between the initial and final states 
of incoming and outgoing radiation. In most works, the Hawking effect could be recovered 
at leading order for some ranges of the black hole’s temperature and the particle’s frequency 
under some suitable assumptions although the mechanisms for recovering it were different 
than in the non-dispersive case. These assumptions and ranges in the models have been 
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briefly reviewed in 18,119]. While most works used the free-fall preferred frame to modify 


the dispersion relation, accelerated frames has been studied in 


111 ]. It has been found that 


the fluxes emitted by a black hole were not significantly affected by the choice of preferred 
frame as long as the acceleration of the frame was not too drastic. However, it has been 
numerically shown that the created particle flux dropped significantly if the acceleration 
became large. Thus, it was conjectured that there might be no Hawking radiation for a 
static frame. 

There are various methods for deriving the Hawking radiation and calculating its temper¬ 
ature. Among them is a semiclassical method of modeling Hawking radiation as a tunneling 
process. This method was first proposed by Kraus and Wilczek 20, 2l|, which is known as 
the null geodesic method. They employed the dynamical geometry approach to calculate 
the imaginary part of the action for the tunneling process of s-wave emission across the 
horizon and related it to the Boltzmann factor for the emission at the Hawking tempera¬ 
ture. Later, the tunneling behaviors of particles were investigated using the Hamilton-Jacobi 
method 22-f24]. In the Hamilton-Jacobi method, one ignores the self-gravitation of emitted 
particles and assumes that its action satisfies the relativistic Hamilton-Jacobi equation. The 
tunneling probability for the classically forbidden trajectory from inside to outside the hori¬ 
zon is obtained by using the Hamilton-Jacobi equation to calculate the imaginary part of 
the action for the tunneling process. Using the null geodesic method and Hamilton-Jacobi 
method, much fruit has been achieved 25436|. Furthermore, the effects of quantum gravity 


on the Hawking radiation have been discussed in the Hamilton-Jacobi method. In fact, the 
minimal length deformed Hamilton-Jacobi equation for fermions in curved spacetime have 


been introduced and the modified Hawking temperatures have been derived 37^42]. These 


motivate us to use the Hamilton-Jacobi method to study the dispersive field theory models. 
Note that a tunneling model has been introduced in dispersive models of analogue gravity 
in [it]]. The Hamilton-Jacobi equations were imported to curved spacetime using the static 
preferred frame in 371-142], which leads us first to considering the static preferred frame for 


the dispersive field theory models. The models with free-fall preferred frame wi 


37 


1 be inves- 


4lj will be 


tigated in 43]. Comparisons between the results in our paper and those in 
given at the end of the section [Til 

The remainder of our paper is organized as follows. In section HH the deformed Hamilton- 
Jacobi equations are derived for the dispersive models. In section ITVl we solve the deformed 
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Hamilton-Jacobi equations to obtain tunneling rates for massive and charged particles to 
O (mp 2 ) and massless and neutral particles to all orders. Thermodynamics of radiations 
near the horizon is discussed in section m The thermal entropy of a massless scalar held 
near the horizon is computed in section [V] using the brick wall model. In section I VII we 
calculate luminosities of a 4D spherically symmetric black hole with the mass M rn p and a 
2 D one. Section IVlII is devoted to our discussion and conclusion, where the limitations of our 
calculations are discussed. Effective held theories incorporating the MDR are constructed 
in the appendix to obtain the deformed Hamilton-Jacobi equations. Throughout the paper 
we take Geometrized units c — G — 1, where the Planck constant H is square of the Planck 
Mass m p . 


II. DEFORMED HAMILTON-JACOBI EQUATION 


In most cases, the modified dispersion relation (MDR) could take the form of 

2 2 tt ( E m \ 

P = ™ P H — i — : 
p \m p rripj 


( 2 ) 


where m p is Planck mass and H (x, y) = x 2 — y 2 for the unmodihed dispersion relation. 
Taylor expanding the right-hand side of eqn. (J2]) for E , m rn p gives 


2 GA E' l m l 
P = K j i+j-2- 


i,j =0 


m v 


( 3 ) 


where h t j is the coefficient of x l yi in the Taylor series of H (x,y) evaluated at (0,0). Since 
eqn. d 3 J) has to become eqn. ([!]) when m p —> oo, we find 


h o,o — ho,i — h\,o — hi t i — 0 and h. 2,0 — ^ 0,2 — 1 - 


( 4 ) 


After some manipulations, eqn. (|3]) can be put in the form of 


p = a 


— ) E 2 - f3 ( — ) m 2 + 7 f —) mE + V 
m„ / \ m r , / \ m pJ ^ 


\ m j 


CL (^)E 


n> 3 


m. 


n —2 


( 5 ) 


where a (0) = (3 (0) = 1 and 7 (0) = 0. If the modifications to the dispersion relation are 
suppressed by some the scale of Lorentz violation A, the naturalness in effective held theories 
would imply that C n ~ (^-) n . For A -C m p , C n could become much large. To include a 
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broader class, we consider a static black hole in the possible presence of electromagnetic 
potential A M with the line element 

ds 2 = f (r) dt 2 — dr 2 — C (r 2 ) h a j } (x) dx a dx b , (6) 

where / (r) has a simple zero at r = with f ( 77 ,) being finite and nonzero. The vanishing 
of / (r) at point r = 77 , indicates the presence of an event horizon. We also assume that the 
vector potential A M is given by 

Ap = A t (r) 60, ( 7 ) 

which is true for charged static black holes in most cases. 

The MDR breaks the Lorentz invariance in flat spacetime. Thus, one needs to pick up a 
preferred frame to determine the form of the MDR. The energy and momentum in eqn. ([5]) 
are defined with respect to the preferred frame, where can be described by the unit vector 
u M tangent to the observers’ world lines. Explicitly, we have 

E = p^, 

p 2 = E 2 - p'jf, ( 8 ) 

where p^ is the energy-momentum vector and E and p are the energy and the norm of 
the momentum measured in the preferred reference frame, respectively. When introducing 
the MDR into curved spacetime, we use the vector held # (i„). To obtain the deformed 
Hamilton-Jacobi equation incorporating the MDR, it is necessary to specify the profile of 
the preferred frame in the black hole spacetime. One of natural frames is a static frame 
hovering above the black hole. For such a frame, the vector field u M (x u ) is 

«"W=(v7 5 ,o) = ^- 7 ==,oj. (9) 

Plugging eqn. (J5J) into eqns. (jHJ), one finds that the energy and the magnitude of the 
momentum becomes 

p = Pt 

vm' 

p 2 = f (r) p 2 r + h c ^ PaPb- (10) 

It can be shown that, if the classical action / is a solution of the Hamilton-Jacobi equation, 
then the transformation equations give 

Pm = -dpi, (11) 
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where — appears since pp = (. E , —p) in our metric signature. Furthermore, since dt is a 
Killing vector of the background spacetime, (cA)^ p^ = Pt is a constant. In fact, p t is the 
conserved energy of the particle and we define uj = p t = —d t I, which means we can separate 
t from other variables when solving for /. Relating / to p M via eqn. (ITT]) and putting 
eqns. (ITU into eqn. (J5J) give the deformed Hamilton-Jacobi equation. In the appendix, the 
deformed Hamilton-Jacobi equation is also derived in a more rigorous way, specifically in 
the language of the effective field theory. We show there that if a scalar/fermion obeys the 
MDR given in eqn. ([5]) in flat spacetime, the deformed scalar/fermionic Hamilton-Jacobi 
equation with respect to the preferred static frame in the black hole background spacetime 
can be both written as 


X 2 = a 


m 


m. 


-b( 


T — ft 


m 




m. 


TTl + 7 


m 


771, 


mT + ^ 


a 


n> 3 


m, 


n —2 


( 12 ) 


where we define 


T 


d t I + qA t 2 

v/TM ’ 


f{r) ( d r lf + 


h ab (x) d a Id b I 
C (r 2 ) 


is the black hole’s electromagnetic potential and q is the particle’s charge. 


(13) 


III. TUNNELING RATE 


In this section, we use the Hamilton-Jacobi method to investigate the particles’ tunneling 
across the event horizon r = r^ of the metric ([6]) by solving eqn. ca. Taking into account 
dtI = — oo, we can employ the following ansatz for the action / 


/ = -ojt + W (r) + 0 (x ), 


(14) 


where uj is the particle’s energy. Plugging the ansatz into eqn. we have 

cu - qA t (r) 


T = 


VJV) 


X 1 = f (r) [d r W(r)} z + 


2 h ab (x) d a Q (x) db<3 (x) 


C (r 2 ) 


(15) 


The method of separation of variables gives the differential equation for 0 (x) 


h ab (x) d a O (x) d b O (x) = A, 


(16) 
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where is A is a constant and determined by h ab (x). Thus, one has 


= / (r) [&W (r )] 2 + —U (17) 

and eqn. (TT 2 l) becomes an ordinary differential equation for W (r). In this section, we solve 
eqn. (fl2l) for d r W (r), calculate its residue at r = r^ and hnd the imaginary part of / which 
gives the tunneling rate T across the event horizon. We will calculate Im W for two cases, a 
massive and charged particle to O (m“ 2 ) and a neutral and massless particle to all orders. 


A. Massive and Charged Particle to O (rn p 2 ) 

Consider a particle with the mass m and the charge q. Solving eqn. (TT 2 l) for p r = d r W (r) 
gives 


Pr = 


±1 


a 


m \ lj (r) 


a/TF) V \ m pJ f( r ) 


~P 


771 


m r 


A 


m — 


C (r 2 ) 


+ my 


m \ u r 


m r 




( 


1 + 




a 1 m 

■p 


Cu n (r) 


\ 


^(r) 


m p~ 2 f (r) 


, ( 18 ) 




where +/— denotes the outgoing/ingoing solutions and lu (r) = co — qA t (r). Here, we have 
a pole at r = ry. Using the residue theory for the semi circleH, we get 


Im W± (r) = ± 71 (i + A qm + O (m p 3 )) , 


2 K 


(19) 


where we define k — f (rh) /2 and 


x Csmy , 1 

t—^am . ~~T 


• !///,/> 2 32/// 2 /.'n 1 


(24C 4 a - 6 Of) J (r h ) uj (r h ) a 2 + Km 2 7 2 i^-i^A 


—koi (6C 2 — 8C 4 a) ( ^ 2 + /3m 2 ^j + a 2 (C ; 2 — 4C 4 a) — ^ ^ ^ 


C(r 2 ) 

The argument is suppressed for a @ 


K 


( 20 ) 


, 7 


and 


in eqns. (1T9|) and 


[1] This procedure will be discussed in detail later in this section. 
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B. Massless and Neutral Particle to All Orders 


We now work with a particle with m = 0 and q — 0. Solving eqn. (fT2]l for p r gives 


Pr = 


±1 


or 


A 


V7MV/W C(r 2 ) 


1 + 


-E c ” 


00 


'771, 


n —2 


/ M’ 


( 21 ) 


/(r) C(r 2 ) n>3 P 

where = C n (0) and we use a (0) = 1. To get the residue of p± at r = r/j, we first define 
a few coefficients C“, , and r] k as follows 

(i+n“ = E c " 


ya n 


71>0 


E c - 

\n =0 


+3^ 


J2C m , n x n , 


n =0 
k 


Pi = 


v—m+i, 


53(-l)'C&<7, (22) 

m =0 

where m is a non-negative integer and C“ are generalized binomial coefficients with C“ = 

n 

a ~fc +1 • Therefore, one has 

—m+| / oo 


fc=l 


= 


±E 


oo 




™ 0 ) / ( r ) 


m=0 P 

oo oo 


f(r)X 

C (r 2 ) oh 


q+3 

^ / f M 


E 


= ± 


co 


EE 


\/ . .k—2l 

A rh uj 


(r 2 ) mp f\~ l ( r ) 


± 


(U 


/ M 




EE 

i=0 fc=0 

A 


A z 


p 2fe+2Z o; 2fc 


C Y/ (r 2 ) m 2fc+2Z f k (r) 


«=o 


^ ^(r 2 )™ 2 


l oo 


E'f 

k=0 


2k+2l ^ k (C\rl) 1 


m 2fc (r 2 ) f k+1 (r) / 


(23) 


where we only keep terms contributing to the residue and set k —> k + 21 in the third line. 

C l (r 2 ) 

Furthermore, we denote the residue of c t^ /fc+i( r ) at r = rh by 

Res ( C ' < r l 1 A = if 

\ C l (r 2 ) f k+1 (r) ’ r / 2k 
Using the residue theory for the semi circle, one has 


urn 


Im W± ( r) = ±— (1 + A), 

ZiKj 


where 


A = E 


l+k> 1 


C(r 2 )m 2 


2 fc+2i /-l 

Vi Si 


00 


2k 


k 


(24) 

(25) 

(26) 


Note that t/{[ — Co — 1 
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C. Calculating A 


It is easy to see that A depends on h ab (x) . Here we consider two kinds of black holes, 4D 
cylindrically and spherically symmetric black holes. For a 4D cylindrically symmetric black 
hole, we have 

h ab (x ) dx a dx b = d6 2 + a 2 dz 2 , (27) 

where — oo < z < oo, 0 < 6 < 2tt and a is some constant. Since dg and d z are the Killing 
holds of the background spacetime, we can separate the variables and consider a solution 
for eqn. (1T6|) of the form 

@ = Jgd + J z z, (28) 


where Jg and J z are constant and Jg is the angular momentum along 2 -axis. The periodicity 
of 0 gives Jg = nh with n 6 Z. Thus, one finds 


A 


J; 


Jz 
H—o ’ 
a z 


(29) 


where Jg = nh with n e Z . 

For a 4D spherically symmetric black hole, we have 


h a b (x) dx a dx b = d6 2 + sin 2 6deft 2 , 


(30) 


where 0 < 9 < n and 0 < (j) < 2 ti. d$ is the Killing vector so we consider a solution for eqn. 
(1T6]) of the form 

0 = Y {6) + J^>, (31) 


where J^ is the angular momentum along 2 -axis. The periodicity of <fi gives J^ = mh with 
m G Z . Since the magnitude of the angular momentum of the particle L can be expressed 
in terms of pg = dgY (0) and J^, 



(32) 


eqn. (fl 6 l) gives A = L 2 . Putting eqn. (j3Tj) into eqn. (TT 6 T) . one gets 


pg = \ A - 


m 2 h 2 
sin 2 6' 


On the other hand, the Sommerfeld quantization for pg gives 



2t r 



ft, 


(33) 


(34) 
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where n is a non-negative integer. The integral in eqn. (1341) is calculated in the classically 
allowed region where pg is real, which requires that A > m 2 h 2 . Integrating the quantization 
integral, one finds that eqn. (15T1) becomes 



(35) 


Solving eqn. (I35|) for A gives the WKB leading quantization of the angular momentum 



(36) 


where l — n + \m\ = 0,1, • • • with \m\ < l. 
quantization of the angular momentum L 2 = 
L 2 = (l + \) 2 k 2 is f. 


Note that the difference between the exact 
l (l + 1) h 2 and the WKB leading quantization 


D. Tunneling Rate 


When one ca' 


dilates the quantum tunneling rate from Im W±, there is so called “factor- 


two problem” 44], Thus, one may have a black hole temperature which is twice the expected 


result. One of solutions is proposed by Mitra 45]. Mitra noted that in general, the action / 
could include some complex constant of integration fC. In this way, the imaginary part of / 
becomes 

Im 4 = Im W± + Im K (37) 


+/— denotes the outgoing/ingoing solutions. In the semi-classical method, the absorption 
probability and the emission probability for a black hole are given by 


P emit oc exp 


P abs oc exp 


( 2 A 

1 — — Im I + J = exp 

-? (Im W + + Im/C) 
h 

( 2 A 

1 — — Im /_ J = exp 

— — (Im W- + Im fC) 
n 


(38) 


On the other hand, it is noted that the classical theory of black holes tells us that an 
incoming particle is absorbed with the probability equalling to one. Thus, one can choose /C 
to impose the classical constraint on the absorption probability, which is Im/C = — ImW_. 
So eqn. (|38|) gives that the probability of a particle tunneling from inside to outside the 
horizon is 


P ern it oc exp 



(Im W + 


Im ITT) 


(39) 
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Another way to circumvent this problem is considering both the contributions from spatial 


and temporal parts of the action to the tunneling ra 


Spatial Contribution: As pointed out in 44], 


46 


es. 


47], Im W± = Im f pfdr used in eqns. 


were not invariant under canonical transformations. Instead, one should take the closed 
contour integral j)p r dr = f p+dr — f p~dr, an invariance under canonical transformations, 

for the tunneling rates P M/ai , oc e X p (±± Imi&dr) ■ 

Temporal Contribution: As shown in 4T-f49j], the temporal part contribution came 


from the ’’rotation” which connects the interior region and the exterior region of the black 
hole. It was found in 49] that the direction in which the horizon was crossed did not 
affect the sign of the temporal contribution. However, the sign of the spatial contribution 
changed when the direction was reversed. Thus, the temporal contributions to P em it/abs 
were the same. When the horizon was crossed once, the action / got a contribution of 
Im (a;At) = and for a round trip, the total contribution was Im (a;At) = —. 

Taking into account the spatial and temporal contributions, one has for the absorption 
probability 


P abs oc exp 
and for the emission probability 

Pemit oc exp 


— Im (pp r dr + Im (a;At) + Im /C 


1 

h 




Im ® p r dr + Im (a;At) + Im /C 


(40) 


(41) 


where /C is a constant of integration. Imposing the classical constraint on the absorption 
probability, one gets 

Im /C = Im ® p r dr — Im (a;At), (42) 




and 


Pemit OC 


— — ( Im (pp r dr 




(Im W + — Im W_) 

h 


(43) 


= exp 

Both approaches give the same expression for P em it■ There is a Boltzmann factor in 
Pemit. with an effective temperature. Using eqns. (Q~9j) and (l25|) , we find that the effective 
temperature for a massive and charged particle is 

To 


^ V® (1 + ^qm) 
and that for a massless and neutral particle is 

T 0 


+ O (m p 3 ) , 


Teff = 


1 +A’ 


(44) 


(45) 
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where we define T 0 = and take ks — 1. 


E. Discussion 


In this section, we suppose that outgoing particles tunnel from r\ < to r 2 > while 
ingoing particles from r 2 to r \. To to obtain the imaginary part of / for the tunneling 
process, we have to give an prescription for evaluating the integrals of Im W± = J^ 2 pf dr. 
Following the Feynman’s ie-prescription 50l|. we take the contour of the integral to be an 
infinitesimal semicircle below the pole at r = ry for outgoing particles. Thus, the integral 
becomes 

i rrh~£ r rr 2 

p+dr — I p+dr + / p+dr+ \ p+dr, (46) 


rr 1 r r h~ £ 

C. 1 1 1 1 w/i 1 1 

1 ri J ri d Cb,£ dr h +£ 


where we denote the semicircle centered at r = with the radius of i? going from + 
— i? to rh — R/rh + R in the upper/lower half complex plane by Cjj/b,r ■ Since the 
contributions from the ranges (r 1; — e) and (r^ + e, r 2 ) are real, the imaginary part of W + 


is 

Im W + = Im / dr. 

Jc B ,e 

Similarly, one has for ingoing particles 

Im ITT = Im / p~dr. 

Jc U:£ 


(47) 


(48) 


To get Im W ±, we expand pp in powers of 


Pr=Yl 


\J f{r)m p 

u> n+1 pi ( r ) 


(49) 


n > 0 / = +1 (r)mf 

where p± (r) are some analytic functions of r around r = rv For the non-dispersive case, 
only the first term in eqn. (1491) appears. Thus, we can Laurent expand / (r) with respect 
to r at r = r h to evaluate f Cg/u pfdr as e —* 0. However, these expansions for pf in the 


dispersive models look suspicious on Cjj/b,£ as e —* 0. In fact, 


’/(» 


can become larger 


than 1 if r is close enough to rp. Thus, we can not trust the expansions for pf any more on 
Cu/B, e - Nevertheless, we can assume that the singularity structure of p r in the dispersive 
models is the same as that in the non-dispersive case except the order of the pole at r = 77 ^ 
which means the MDR effects do not introduce branch cuts or new poles for p r in the upper 
or lower half complex plane. Note that one may need a complete theory of quantum gravity 
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to justify the assumption. Now consider the semicircles Cu/b,r with large enough R, which 
lies in the region where the expansion for p r can be trusted. Under the assumption, there 
are no poles inside the area enclosed by (R, r h — e), Cu/B,ei ( r h + e, R), and Cu/b,r- Thus, 
we have 


rrh-R 


Im W± = Im 


pfdr + Im 


pfdr + Im 


rr 2 


pfdr 


'r i 


Ic 


B/U,R 
n +1 


f o; n+1 f 

= Im / ptdr = > -Im / 

/ r ' 77l n In 

J ^B/U,R n> 0 P J( -'B/U,R 


r h +R 
r) 


Pn 


c B/UM f^ +1 (0 


dr, 


(50) 


where contributions from the ranges (r^r^ — R) and (r^ + R, r 2 ) are discarded since they 
are always real. If the radii of the Laurent series of y expanded at r = rh are larger 


-l 


than R, we can Laurent expand ^ on Cb/u,r an d only the coefficients a_i of (r — r h ) 

terms contribute to the imaginary part of the integrals. This justifies the procedure to 

obtain to obtain eqns. (TT9j) and fj25jl . Since the expansions for pf can be trusted on Cjj/b,r , 

2 

< 1 on Cu/b,Ri which gives R > - i U I . Usually, one has that the radii of 


\J uRn 




the Laurent series of 


Pn(U 


~ K 


-1 


and hence k 1 > R. For R to exist, one has k 1 > 


^ K,rrip 


which leads to u < m p . Note that if Lorentz violating scale A is much smaller than m p , one 


instead has that 


'/(*■) a 


< 1 on C u/b ,r and u < A. 


Various theories of quantum gravity, such as string theory^ loop quantum gravity and 
quantum geometry, predict the existence of a minimal length 51-153]. The generalized un¬ 
certainty principle(GUP) 54] is a simply way to realize this minimal length. To incorporate 
the Klcin-Gordon/Dirac equation with the GUP, one usually considers the quantization in 
position representation. In position representation, the operators k = —iS7 and oj = id t are 
introduced 55]. One then can express the energy and momentum operators as functions of k 
and oj and obtain the deformed Klein-Gordon/Dirac equations in flat spacetime. Inserting 
the ansatz (p = exp ( iEt — ip - x) in the Klcin-Gordon/Dirac equation gives the dispersion 


relation for E and p in flat spacetime. In 37 4jJ, the deformed Dirac equation was general 


ized to curved spacetime. The modified Hawking temperatures of various black holes were 
then derived via the Hamilton-Jacobi method. It turns out that the way of generalizing 
the Dirac equation to curved spacetime in 37-1411] is the same as that in the appendix of 
the paper where we choose the static frame as the preferred one. Thus, we can use the 
dispersion relation for E and p obtained in flat spacetime and eqn. (1201) to reproduce the 


modified Hawking temperatures of black holes with the metric (J6|) obtained in 


37 


4l|. In 
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fact, the dispersion relation in flat spacetime in 37-41[ is given by 


p 2 « E 2 - m 2 + 2 PE a + £> (/3 2 ) , 


( 51 ) 


which by comparing to eqn. (J5]) gives 


a = 1, f3 = 1 ,7 = 0, C 3 = 0, C/t = 2 , and rn p = —^=. 


1 

VP' 


(52) 


Thus, eqn. 


becomes 

24a;' (r h ) Co (r h ) 


A qm = P 


1 

+ - 


A 


+ rri 


f" (r h ) Co 2 ( r h ) 


(53) 


16k ' 2 (r 2 J ' J 4k 2 

For example, A qm for a particle with the angular momentum l — 0 in a Schwarzschild black 
hole with the mass M[37l. l4l| is 


= C ( m 2 + 4w 2 ) 


(54) 


Aq m for a neutral particle with t 


the mass M and the charge Q 38 


le angular momentum l — 0 in a Reissner-Nordstrom with 


4lj is 


Aqm P 


m 

T + 


2 r i (r i — 2 r 


lo 2 


(55) 


(r + — r _) 2 

where r± — M± a/ M 2 — Q 2 . ffere we only consider a neutral particle to make a comparison 
since the electromagnetic hied was included in 38, [41] in a different way. 

Following the argument proposed in 56], the authors in 57] obtained modified relations 
between the mass of a 4D Schwarzschild black hole and its entropy and temperature. The 
argument connecting a MDR and some modifications of the entropy of black holes is for¬ 
mulated in a scheme of analysis first introduced by Bekenstein 
temperature of the black hole for the MDR in eqn. (J5l) with m = 0 was given by 


58]. In fact, the modified 


T = Tn 


1 - 


m p C 3 


+ 


m; 


4 M AM 2 \ 8 


7 C 2 C d 


+ 0 1 p 


m: 

M 3 


(56) 


where M is the mass of the black hole and T 0 = Aka. On the other hand, we can use 


eqn. (l20l) to estimate the temperature of the black hole. For a massless particle in a 4D 
Schwarzschild black hole, eqn. (l 20 l) gives 


A = 


32m 2 L 


A 


(4C, - 3 Cl) — 2 + 8 (4C, - Cl) 'J 


+ O (m - 3 ), 


(57) 
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where A is the magnitude of the angular momentum of the particle. The event horizon of 
the Schwarzschild black hole is = 2 M. Near the horizon of the the black hole, one has 
A r ''“' ( pr h f (cnr/i) 2 . Thus, one can rewrite A 


A 


cu 


16 m^ 


( 2 OC 4 - 7C 3 2 ) + O (m“ 3 ) 


(58) 


As reported in |57|, the authors obtained the relation u > ^ + O between the energy 
of a particle and its position uncertainty for a MDR. Near the horizon of the Schwarzschild 
black hole, the position uncertainty of a particle is of the order of the Schwarzschild radius 
of the black holejssj Sx ~ = 2 M. Thus, one finds for T 


T ~ Tn 


1 + 


mi 


64 M 2 


(7Cl - 20C 4 ) + O 


M 3 


(59) 


Note that the term — in eqn. (I55j) could imply that the leading correction to the entropy 
of the black hole should have y/area dependence. Since in Loop Quantum Gravity such a 
y/area contribution to black-hole entropy has already been excluded, C 3 should vanish 56]. 


However, it has been shown in 


57] that this might be evaded by combining MDR and GUP. 


Our result suggests that a nonvanishing C 3 in the MDR does not lead to the presence of 
y/area contribution to the black hole entropy. 


IV. THERMODYNAMICS OF RADIATIONS 

For particles emitted in a wave mode labelled by energy u and A plus some other quantum 
numbers J t if needed, we find that 

(Probability for a black hole to emit a particle in this mode) 

= exp ( — —— J x (Probability for a black hole to absorb a particle in the same mode), 

V T effJ 

where T e // is given by eqns. (j44l) or (1451) . The above relation for usual dispersion relation 
was obtained by Hartle and Hawking|59l] using semiclassical analysis. If the black hole is in 
equilibrium, the rate of emission particles by the black hole must exactly equal the rate of 
absorption. Neglecting back-reaction, detailed balance condition requires that the ratio of 
the probability of having N particles in a particular mode with oj, A and J* to the probability 
of having N —1 particles in the same mode is exp (^ — T 77 ) • Thus, we find that the probability 
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of having N particles Pn (a/, A,Jj) in the mode is given by 


P N ( u , A, Ji) = C uXJi exp 


Nu 


T, 


eff 


(60) 


where C uX j i is a normalizing constant. C U) \ y j i is determined by the normalized condition 


N a 


I> (a/, A, Jj) = 1 where iVoo = oo for bosons and iV^ = 1 for fermions. Thus, the 


w=o 


probability Pn (u, A, Jj) is 
Pn (n-> A, Ji) 


1 - (—l) e exp 


U 


T f 


eff 


-i 1—2e 


exp 


Nu 


T f 


eff 


(61) 


where e = 0 for bosons and e = 1 for fermions. To calculate the average number n u> x,Ji i n 
the mode, we define 


N a 


(v) = exp ( N v - 


N =0 


Nu 


eff 


(62) 


where one has C LJi \,j t = \ 3 . (0). So we find 


N a 


^ ] NPn (u, A, Ji 


N =0 


(j) 

(j) 


1/2=0 


exp ) - (-1) 


(63) 


Using eqns. ()6T]) and (163|) . one can rewrite Pn {ui, A, Jj) in terms of n uX j i as 

Pn (u, A, J f ) = [1 + (-1)' , 


(64) 


where N can be any non-negative integer for bosons (e = 0) but is restricted to be 0 or 1 
for fermions (e = 1). The von Neumann entropy for the mode is 

Woe 

Su>,\,Ji ^ ) Pn (^h A, Jj) In Pn (w, A, Jj), 

w=o 

= [n u ,x,ji + (-1)1 In [1 + (—l) e n Ut x,ji\ - \nn uXJi (65) 

Woo 

where we use (u, A, Jj) = riu,x,j r The total entropy of radiation is 

w=o 

S = Va ,j v (66) 

LJ,X ,Ji 

which will be calculated in the brick wall model in section [V] Note that since T e ff only 
depends on u and A, the average number n^x,^ and the entropy s^a.j* are independent of 
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Jj. Thus, we could omit the subscript J* in n ,and from now on. Defining n (x) 

and s ( x ) by 


n ( x ) = 


s (x) = 


exp x — (—l) e 
(—l) e exp x 


exp x 


- 1 )' 


In 


exp x 


expx 


-ir 


+ 


In [exp x — (—l) e ] 


exp x 


■i)' 


(67) 


we can write n U) \ and s u ^\ with respect to n ( x ) and s ( x ) 


?Va = n 


Suj,\ 5 


C 0 


T ( 


eff 
C0 


T, 


eff 


( 68 ) 


When integrating over u, we need to specify the upper limit on the energy of the emitted 
particle. One of the limits comes from the requirement that the energy of the particle could 
not exceed the mass of the black hole. Another one is from the effective field theories in the 
appendix. Suppose that the higher dimensional operators in the effective held theories are 
suppressed by some scale of Lorentz violation A. Usually, we can only trust the effective 
theories below A. As decoupling theorem 60] shows, the contributions above A in some 
regularized theory gets absorbed into Wilson coefficients of the effective theories, Cs and 
Bs in the appendix. Consequently, the energy of the particle could not exceed A otherwise 
our effective theories would break down. Note that cu < A has also been obtained in section 
m Thus, the energy of the emitted particle to < cu max = min{M, A}. In the remaining of 
the paper, we would encounter the integrals like 


u l n (u) du or 


u l s 


u 


du , 


(69) 


Jo Jo 

where M max = and i is a non-negative integer. For a black hole with the mass M m p , 
one hnds that w max = ~ 1. For example, k = jjt and Km v 1 in the 

Trip AC KTYlp etivj •* 

Schwarzschild metric. For such case, using n (x) ~ e~ x and s (x) ~ xe~ x for i> 1, one gets 


udn (u) du 


1 

KTUrt 


[KTYlr 


and 


u l s 


(u) du 


1 

Km?) 




\i+l ’ 


(70) 


which can be safely neglected for nm p 1 and hence we can let u max = oo in eqn. flfFUl . 
Therefore, in section [V] we let u max = oo for integrals of the form in eqn. (1691) since we are 
only interested in the divergent part of the entanglement entropy as nm p —> 0. 
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Most works of the dispersive models focus on calculating the modifications of the asymp¬ 
totic spectrum. There are two types of contributions to the asymptotic spectrum measured 
at infinity: particles generated at the black hole horizon and scattering off the background. 
The first contribution depends on derivatives of / (r) at r = whereas the second one 
depends on / (r) for r > r^. As shown in section IVll the the spectrum is proportional to 
n U) \. Since the tunneling probability across the horizon is calculated in the Hamilton-Jacobi 
method, only the first contribution to the asymptotic spectrum is computed in our paper. 
The dispersive models have been studied in 2D spacetime for most works. The higher or¬ 
der terms in the MDR violate conformal invariance of 2D spacetime, hence there is some 
scattering. Our calculations show that the spectrum of created particles near the horizon 
is close to a perfect thermal spectrum in the dispersive models with the static preferred 
frame, where it has been suggested that the asymptotic spectrum could significantly differ 
from the thermal one 11], Thus, one might need the scattering off the background to ex¬ 


plain the difference between the spectrum near the horizon and the asymptotic one. The 

near horizon properties of the dispersive modes has been considered in 16]. It was found 

that the mode mixing responsible for the Hawking effect occurred within a region of size 
— 1 . 

db r ~ rripK~ 3 across the horizon and the standard Hawking spectrum has been recovered as 
long as db r -C k -1 . 

In section [V] we calculate the thermal entropy near the horizon, where the effects of 
scattering off the background hardly play a role and hence are neglected. However, the black 
holes’ luminosities at infinity are computed in section m where the effects of scattering off 
the background should be included. In fact, we use the geometric optics approximation to 
estimate these effects instead of solving the equations of motion of the fields. Under the 
geometric optics approximation, we find that the asymptotic spectrum is still close to a 
thermal one. Nevertheless, [ll| might suggest that such approximation may not be robust 
enough and hence the asymptotic spectrum could be significantly changed. 


V. ENTROPY IN BRICK WALL MODEL 


In 1985 t’ Hooft 61] proposed the brick wall model to calculate the entropy of a thermal 
gas of Hawking particles propagating just outside the black hole horizon. The entropy is 
calculated by methods of the WKB approximation. However, when it comes to calculate 
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the density of states of emitted particles, t’ Hooft found that they became infinite as one 
got closer to the horizon. To make the entropy finite, he introduced a brick wall cut-off near 
the horizon such that the boundary condition 

$ (x) — 0 at r = r h + r e , (71) 


where <J> is the radiation’s field. Moreover, another cut-off at a large distance from the 
horizon L>r/, was introduced to eliminate infrared divergences. 

For simplicity, we consider in this section a 4D spherically symmetric black hole with the 
metric 

ds 2 = f (r) dt 2 ——dr 2 — C (r 2 ) dfl. (72) 

/ (H 

For such a black hole, the quantum numbers needed to specify a wave mode of radiation are 
the energy u, the angular momentum l, the magnetic quantum number m. We also assume 
that the radiated particles are massless and neutral. Thus, the MDR, eqn. (J2D, becomes 

p2=m K^°)= m Kw)' (73) 

where we define H (x) = H ( x , 0) and H (x) ~ x 2 for x <C 1. As shown in section m 
the deformed Hamilton-Jacobi equation incorporating eqn. (1741) for a massless and neutral 
particle in the 4D spherically symmetric black hole is given by 

* 2 = < H (£) • < 74 > 


where 


T = 

Then, we get from eqn. (174)1 


cu 






/ (r )pi + 


{i + \) 2 n 2 

C (r 2 ) 


(75) 


Pr = 


f(r) 


m 2 H 


UJ 


m. 




Dehne the radial wave number k (r, l,u) by 


(/ + |) 2 h 2 

C (r 2 ) 


(76) 


k(r,l,u) — \p r \ , (77) 

as long as p 2 > 0, and k (r, /, to) = 0 otherwise. Taking two Dirichlet conditions at r = rh + r E 
and r = L into account, one hnds that the number of one-particle states not exceeding c o 
with fixed value of the angular momentum l is given by 

n(cu,Z) = — [ k(r,l,oj)dr. (78) 

7r "' Jr h +r e 
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Thus, we obtain for the total entropy of radiation 


S=^2 s Ut i = I (21 + 1 ) dl J du 


( 2 1 + 1 ) dl [ du [ L dr dk ^ l,UJ K u , 
J Jr h +r e A "' 


l 

7T h _ 

c (rp m 2 p r L dr 

Ir h +r e f i r ) 

c( r r ' 


dn (u, l ) 
du 




rH 


2nh? 

m P Vm) C (rl) ^ 


duH’ 


du 

u 


m 


H 


u 


m r 


— z 


C{rp 

C (r 2 ) 




(79) 


where we define a dimensionless parameter z = . Using A = {l + \) 2 K 2 and 0 = 






c(, 


h) m l 


, one rewrites eqn. 


as 


A = 


OO CO 

EE 

a =0 k =0 


2k+2a /-a 
'la 


U 


2k 


m, 


2k 


- 1. 


(80) 


Defining the coefficients k by 


one has for 


CO CO 

<,k'+2i'fV j 

Z'=0k'=0 


Ck'Z l 'x 2k ' - 1 


= EE«U*“ 

fc=0t=0 


l s 


u 


T, 


eff 


= .,'" (1 _ +A) ) 


(n) 


n! 


n=0 
oo co co r 


E 

n=0 
oo 

EEE 

a=0 n=0 /c=0 
oo 

(u), 

a =0 


T n 


u \ u n A r 



T 0 ) Tq 


s W (it) iT 


rt! 


zu 



(81) 


(82) 


(83) 


where we define u — and 


k =o 


T 0 u 


rn r 


2k k+a 

v 

s (it) U 


n\ 

n =0 


0 and il<n-k,k = 


f n 

s a.b 


(84) 


Note that l = 0 except = 1 ,f" fc<n _ z = 0 and i) <n _ k j. = 0. Putting eqn. (JH2J into eqn. 
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(179]) and integrating eqn. d79j) over z gives 

5 C(ri)mjT 0 jk ^Fa! 


2irh 3 


r (a + I) y 0 


[ L dr C a+1 (r 2 ) a +i 
K h+ r e f (r) C^ (rl) 


@ a (n )du 
T 0 u 




H' 


T 0 u 


m. 


•v / 77 


(85) 


To calculate S', the variable r may be changed by introducing x = 


Tqu 


becomes 

5 C (rl) ml ^ v^Fa! 


oy/JiF) 


Then, eqn. 


0 a (u ) du / x^H^ (x) 77' (x) G, 


4^2 ^ T (a + |) Jo 


x 2, m^ 


dx, ( 86 ) 


where we dehne x £ = - G= 

m PV/( r h+ r e) 


,5 = 


Tqu 1 


and 


Ga (y) = 


l P y/f(L) 

2kC° +1 [r 1 [yf] 


(87) 


^ a+1 ( rl ) f [/- 1 ( 2 /)]' 

Now the brick walls are at x = x £ and x = 6. Note that x £ —> oo when r £ —> 0 and the 
horizon is at x £ = oo. Since G a (0) = 1, we can Taylor expand G„ ( y ) at y — 0 


Ga ( V ) = 


( 88 ) 


k =0 


where we hnd the first two coefficients of the series expansion are /q = 1 and /“ = 


r h C'(rl) 

KC ( r h ) 


(a + 1 ) - 


C(rpf"(r h ) 
4 K.r h C'(rl) 


Substituting eqn. 


into eqn. 


gives us 


S 


C (r 2 h ) y- / m p K\ ™ ™ y/FFal/g 
An 2 ml^\2n) ^ T (a + f) 


u 2k Q a (u) du 



77 a+ 2 (x) 77'(x) x 2/c 1 dx, 


(89) 

where we use T 0 = and h = m 2 . The entropy receives two contributions, one from the 
horizon and the other from the vacuum surrounding the system at large distances. The 
second one is irrelevant for our purposes and henceforth discarded. 

For the usual scenario with 77 (x) = x 2 , the integrals over x in eqn. (1891) become divergent 
for a = k = 0 as one approaches the horizon with x £ —> oo. This divergence leads to the 
introduction of the wall near the horizon by t’ Hooft. However, the x-integrals could be 
finite as x £ —> oo for some MDRs. In fact, there are two kinds of MDRs for the integrals to 
be finite. For the first kind of these MDRs, the high energy contributions are suppressed. 
For example, the “all-order MDR” of form 


P 


2 


2 m 2 exp 



(90) 
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was given in the ^-Minkowski noncommutative spacetime in 62], For such a MDR, one has 

H (x) = 2 exp (—x) [cosh (x) — 1], 

H l+1 * (x) H' (x) x~ 2k ~~ l 


exp (—x) x 


—2k—l 


(91) 


as x —» oo, 

which guaranties the convergence of the ^-integrals as x £ —> oo. Another example is inspired 
by the all order generalized uncertainty relation considered in 
as 


dp 

~dE 


= exp — 


E 2 


m. 


63]. The MDR can be written 

(92) 


which gives 


H(x) = 


H l+ 3 (x) H' (x) x-™- 1 < (e~ x 


dx 


21+3 


x as x oo. 


(93) 


Thus, the x-integrals stay finite as x —> cxd. Moreover, it was found in 63] that the entropy 
kept hnite when the wall approached the horizon and hence the wall in the brick wall model 
located just outside the horizon could be avoided. For the second kind of the MDRs, the 
energy E in the MDRs has a maximum value and hence x e can not go to the infinity. For 


example, Corley and Jacobson 


9] proposed 


E = „ Ip 2 


/r 


4 m?p 


(94) 


which gives 


H l+l * (x) H' (x) x 


H (x) = 2 (l — y/l — x) for 0 < x < 1, 
1 


—2k—l 


as x —> 1. 


(95) 


\/l — x 

So the x-integrals are hnite for the Corley and Jacobson dispersion relation, ft was shown in 


64 } that the entropy was rendered UV hnite for the Corley and Jacobson dispersion relation. 
For the Unruh dispersion relation[7] 


E = m r 


tanh I ?- 


m, 


p / j 


(96) 


we have 


H l+ 2 (x) H' (x) x 


H (x) = [tanh 1 (x n )] n for 0 < x < 1, 
[In (1 — x n )]^ + 2)^ _1 


—2k—l 


l-x r 


as x —» 1. 


(97) 
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We find that the ^-integrals diverge as x —> 1 and a wall near the horizon is needed. However, 


the entropy for the Unrnh dispersion relation was also found finite in 


64|. This might be 


due to different generalization of the Hamilton-Jacobi equation in curved space in 64]. In 
the remaining of section, we will consider two cases, in one of which the x-integrals converge, 
and in the other they diverge. 


A. UV Finite Case 

We here assume that the x-integrals converge as x —» A, where A = oo for the first kind 
of the MDRs in this case and A is the largest x for the second kind. Thus, we can define 

r A 

■'k ~ 


? a — 
C7, - 


H a+ * (x) H' (x) x 2 l dx . 


Since H (x) ~ x 2 for x <C 1, the Taylor expansion of H a+ 2 (x) H' ( x ) is given by 


(98) 


H a+1 i (x) H' (x) = d ) x ' 
j =0 


j+2a+2 


(99) 


where <Jq = 2. Then one gets 


c h = 


'XI 


pXl _ _ 

H a+ 3 (x) H’ (x) x-^dx + / d a jX j+2a+2 x~ 2k - l dx 

J S ■ n 


= C, 


E 


3 =0 

^a^+2(a-fc)+2 

9 {k - a - 1) d a 2k _ 2a _ 2 In 5, 


j + 2 (a — k) + 2 


( 100 ) 


j=0,j^2k—2a—2 ' 

where 6 (x) is the Heaviside step function, 0 < Xi < A and c\ is a constant independent of 
L. Neglecting terms depending on L, one finds 


c a k -c a k -9(k-a- 1) d 2k _2 a —2 ^ 


m v KU 
2 vr 


( 101 ) 


Plugging eqn. (11011) into eqn. (189|) gives us that the entropy near the horizon can be written 
of form 


s = kMlf; ( Sk+ i k in 


47r 2 m^ 

p k =0 




2k 


( 102 ) 


2tc J \ 2n 

For k = 0, we can choose X\ = 0 in eqn. (1 101 1) since J Q A H a+ ^ (x) H' (x) x~ 1 dx is convergent 
as x —» 0. Hence, one has 


r a — 

On — 


H a+ 2 (x) H' (x) x 1 dx, 


(103) 
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and 


c“ (ril 


2\ a 


7r 


50 5 r(a + §) ./„ 


s ® ( u ) u l du, 


(104) 


where rjf = ^4- — ^ and we use £" fe<n _; = 0, = (yfCo ) 1 , and Co = 1- Since there is no 
In 6 in eqn. (llOip for k — 0, we have Z 0 — 0. For k = 1, only c® contributes to l\ and we hnd 


where we have 


h = -4/i° 


f0 _ ^c" ( r D 
./1 — 


s (w) u 2 du, 

C ( rl) f" (r h ) 


(105) 


kC (r^) [ 4nr h C' (r 2 h ) 
Thus, we obtains for the entropy near horizon 


S 


C (rg) 

47r 2 m2 


Sq- 


Kr h C (r 2 h ) 
4n 4 


c (■ fj ) r (r h ) 

4nr h C' (r?) 


s (n) u 2 du In m p K+Finite terms as nm p —* 0. 

(106) 


B. Perturbative Case 


In this case, a wall near the horizon is needed to regulate the ai-integrals. As above, the 
function H (x) can be presented in the form of Taylor series 


H(x) = 


X>‘, 


(107) 


i=2 


where C2 = 1. One then can have Taylor expansions for H a+ 2 ( x ) H' (x) 


H a+1 * (x) H' (x) = ^ d)x- 
j =0 


j>+2a+2 


(108) 


where c?q = 2. The radial position of the brick wall near the horizon is r = + r e (x = x e ). 

The invariant distance of the wall from the horizon e is defined by 

f rh+re dr f ye 2 dy 


where we define y = 


£ = 


_ nu = VM an d y = Abm. Noting - 

„ 11. ,y 0 cr.r-m.^, 0 


l rh ^TV) Jo f l (f- 1 (y 2 )Y 

. _„ „ =1, one obtains 

x £ m p 0 2k 


(109) 


(•Vs 


£K = 


dy 


f'if-Hv 2 )) 

2 K 


rye ( _°°_ 

l 1+ 5A 

J 0 \ n=l 


_°°_ f „2n 

y- 1 dy = Vc \ 1 + jr )nV ‘ 


n =1 


2n + 1 / ’ 


( 110 ) 
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where we expand r^Tyj in the integral and f n are coefficients of the series. Solving eqn. 
(1 11 01) for y £ gives 

f OO ^ 

l + ^CnM 2n ) , (111) 


Ve = 


n=l 


where ( n are determined by f n . Using x £ = T ° u , one can relate x £ to e by 


Verrip ■ 


OO 


nrcLrf.a i _ 




n =0 


m p x £ 

In 


T 0 u 


^Xo ( £K ) 2 " ~ ln£K > 


( 112 ) 


71=1 


where Xa — 1- Focusing only on the near horizon contributions, we neglect terms involving 
L and use eqn. (llOSp to obtain 


H a+1 * (x) H' (x) x- 2k ~ l dx ~ 


^0^1+2(a—fc)+2 


j=0 ,j ^2k—2a—2 


■ ,, 9 f0(fc-a-l)d“ fc _ 2a _ 2 ln^, 

j + 2 (a — k) + 2 T 0 n 

(113) 

x £ m p 


where for the logarithmic term we have Jg S x 1 dx = In ^ = In | In / (L) ~ In 

Plugging eqns. (11131) and (|84l) into eqn. (l89l) . one hnds for the entropy near the horizon 
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At hrst sight, it seems impossible to single out the most divergent part of eqn. (II 14)1 since 
j + 2 a in the first term of eqn. (1 114 [) can go to infinity. However, the brick wall is put 
at r = r/j + r £ to cut off some unknown quantum physics of gravity. In this sense, the 
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invariant distance of the wall from the horizon e could be given by e ~ m p . Indeed in the’t 
Hooft’s original calculation for Schwarzschild black holes, requiring that the entropy of the 
radiation near the horizon SBrick the Black hole’s Bekenstein-Hawking entropy Sbh also 

K-m p for a scalar held. Thus, we define a such as £ = am p . Replacing £ by 


gives £ 


> 


907T 


arrip in eqn. (11141) . we find for the entropy 
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where we dehne 
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C. Discussion 


For a massless scalar held, we find the entropy near horizon in both cases can be written 


as 
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167r 3 777,p 


+ sl hi rripK + Finite terms as nm p —> 0, 


(117) 


Kr h C'{rl) 

where s L = - ^— L 


1 - 


C{rj)f"{r h ) 


4 Kr h C'{rl) 

scenario without the MDR, the entropy near horizon 
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r j ) is the horizon area. For the 


65l. 166 ] is 
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where we let the proper distance £ = am p . By comparing eqn. (II 171) with eqn. (I118p 
it shows that the leading term of the entropy is affected by the effects of the MDR while 
the subleading logarithmic term is not. On the other hand, the first law of black hole 
thermodynamics dSs = -fr- and eqn. ()59l) lead to the modified entropy of the black hole 
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where A = I 67 xM 2 and n = yO-. For Sb, the leading term is not changed while the subleading 
logarithmic term is due to the MDR. These might suggest that the explanations for statistical 
origin of the black holes’ entropy need more than the entropy of a thermal gas of Hawking 
particles near the horizon. 
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Since the deformed Hamilton-Jacobi equations and the corrections to the Hawking tem¬ 
perature are same for fermions and scalars with the same MDR, one may wonder if eqns. 
(11061) and (11151) also work for fermions. In fact, it has been shown in 67] that the same 
argument in this section held for fermions if an appropriate boundary condition was taken 
instead of the too restrictive Dirichlet boundary condition. 

For a MDR with H ( x ) in the UV finite case, we have shown that a brick wall near the 
horizon is not needed since the entropy is finite as one approaches the horizon. However, 
if one expands H ( x ) as a power series of x and calculates the entropy in the perturbative 
case, it seems that a wall near the horizon is needed to regulate the divergence. How 
can we reconcile the contradiction? As noted in 68], the divergence of the entropy in the 


perturbative case as a —» 0 is more like due to the breaking down of the Taylor series. 
For the typical energy lu ~ To = , one finds that H ( x ) and the MDR corrections to 


the entropy are powers of 


hx, 




\/ f{r)m p 2ir^/ f{r)n 


At the wall at r £ & + 2nrni, we have 


, - ~ Thus, the perturbative case is valid outside the wall at r e = r h + 2 m7 C. 

v / K r e) rn p 

However, the perturbation would break down deep within the wall and the closed form of 


H (x) is needed. 


VI. BLACK HOLE EVAPORATION 


In 


69], Page counted the number of modes per frequency interval with periodic boundary 


conditions in a large container around the black hole and divided it by the time it takes a 
particle to cross the container. He then related the expected number emitted per mode n 
to the average emission rate per frequency interval ^ by 


dn doj 
dt 2 tt fi 


( 120 ) 


for each mode and frequency interval (w, oj + doj). Following the same argument, we find 
that in the MDR case 


dii duj dp r 

dt dp r 2nh 


= n 


dco 

2nTi 


( 121 ) 


where is the radial velocity of the particle and the number of modes between the wavevec- 

OPr 


tor interval ( p r ,p r + dp r ) is Since each particle carries off the energy the total lumi¬ 
nosity is obtained from multiplying % by the energy to and summing up over all energy uj 
















(122) 


and quantum numbers, denoted by i, 



However, some of the radiation emitted by the horizon might not be able to reach the 
asymptotic region. Before the radiation reaches the distant observer, they must pass the 
curved spacetime around the black hole horizon, which plays the role of a potential barrier. 
This effect on L can be described by a greybody factor from the scattering coefficients of the 
black hole. Actually, the greybody factor is given by \T l (cu)| 2 , where Tj (cu) represents the 
transmission coefficient of the black hole barrier which in general can depend on the energy 
uj and quantum numbers i of the particle. Taking the greybody factor into account, we find 
for the total luminosity 

L = Ylf \ T i( U )\ 2 ( 123 ) 

i 

Since the relevant radiation usually have the energy of order hM _1 , where M is the mass of 
the black hole, one should use the wave equations given in the appendix to compute \T t (cj)| 2 
accurately. However, solving the wave equations for \T t (o;)| 2 could be very complicated. On 
the other hand, one can use the geometric optics approximation to estimate |T ? ; (cj)| 2 . In the 
geometric optics approximation, we assume uj^> M and high energy waves will be absorbed 
unless they are aimed away from the black hole. Hence we have \T t (ca)| 2 = 1 for all the 
classically allowed energy c o and quantum numbers i of the particle, while \Ti(E)\ 2 = 0 
otherwise. For the usual dispersion relations, the Stefan’s law for black holes is obtained 
in this approximation. In the remaining of the section, we will discuss evaporations of a 
4D spherically symmetric black hole with the mass M m p and a 2D black hole. For 
simplicity, we assume that the particles are massless and neutral. 


A. 4D Spherically Symmetric Black Hole 

To find the classically allowed values of angular momentum l with fixed value of energy 
uj, we consider eqn. (1761) for a massless particle in a 4D spherically symmetric black hole, 
where we have A = (l + |) 2 H 2 . Since p r is always a real number in the geometric optics 
approximation, one has an upper bound on A 

UJ 

™ P Vf ( r ) 


A <C (r 2 ) m\H 
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Suppose C (r 2 ) m 2 H 


has a minimum at r m ; n and this minimum is denoted by 


P \m p y/f(r) / 

A m ax' If the particles overcome the angular momentum barrier and get absorbed by the 
black hole, one must have A < A max . Thus, the eqn. (| 123ft becomes 
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where g s is the number of polarization, z = an d we use ec l n - (EHD for Defining 
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where is given by eqn. (IHTlh Substituting eqn. (11271) into eqn. (j 125ft and integrating 
eqn. (1125ft over z gives 
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where we let u max = oo for M ^ m p . Since H (x) = x 2 + C n x n 1 we define h ^ by 
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where hg = 1, = aC 3 and = C 4 a + DiItJE. Plugging eqns. (1129ft and (1126ft into eqn. 

(I128p gives 
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where [x] = max {m G Z \ m < x} and we define 

f°° n W ( u) u j+l+3 C a+1 (rT n ) 
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We now use eqn. (I130j> to calculate the luminosity in the Schwarzschild metric to O 


M- 


For the Schwarzschild metric, one has / (r) = 1 — —, = 2 M, C (r 2 ) = r 2 and k = -jjj. 


Taking the derivative of C (r 2 ) m 2 H 


■>\Jfi r ) 


and equating it to zero, one finds 
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(132) 

,p a 7rt“ V'^p// 

For emitting 77 , species of massless scalars and 77/ species of massless spin-1/2 fermions from 
a Schwarzschild black hole into empty space, putting eqns. (j 132ft into eqn. (1130ft gives the 
total luminosity 


L = 


9 ml 


409607rM 2 


T \ 777, 

n s + -n f J + (0.2671, + 0.5077/) C 3 —j 


171? f 171? \ 1 

+ [(0.1577, + 0.2977/) C 2 - (0.2477, + 0.4677/) C 4 ]j^ + 0 (^J j . (133) 

In the geometric optics approximation, the Schwarzschild black hole can be described as a 
black sphere for absorbing particles. The total luminosity are determined by the radius of 
the black sphere R and the temperature of the black hole T. Note that one has R = yf 
and T e ff ~ T 0 (1 — A) where for massless particles 


A = 


32m 2 L 


(4G, “ 3C h 2M? + 8 ( 4C4 - C ») 


(134) 


Consider a sub-luminal case with C 3 = 0 and C 4 > 0, where the total luminosity decreases 
due to the MDR effects. In this case, the MDR effects increase the radius of the black 
sphere while they decrease the temperature of the black hole. The competition between the 
increased radius and the decreased temperature determines whether the luminosity would 
increase or decrease. It appears form eqn. (1133ft that the effects of decreased temperature 
wins the competition. 

B. 2D Black Hole 

Suppose the metric of a 2D black hole is given by 


ds 2 = f (r) dt 2 - ——dr 2 , 

/( r ) 


( 135 ) 
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where / (r) has a simple zero at r = tv Here we consider a neutral and massless scalar 
particle governed by the modified dispersion relation 

Cp 4 


E —p — 


mt 


(136) 


which is the Corley and Jacobson dispersion relation for C > 0 [[9]. Expressing p in terms 
of E gives 


p 2 = E 2 


CE 4 


mt 


-(S' 


For the 2D black hole with the event horizon at r = rh, eqn. (j20l) gives 
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where we use m = 0, Cj (r h ) = uj, A = 0, a = 1, 7 = 0, C 3 = 0, and C 4 = C. For 


u < l u max , the term in eqn. (1138P dominates and hence the terms O 


are neglected 


for simplicity. Define rj = — which becomes 1 for a 2D Schwarzschild black hole with 

f (r) — 1 — —. In this case, we can choose the cutoff of the effective theories A = Q 1A_ 

r 

where 0 < a < 1. Note that |A| < 1 for cu < A. Therefore, the luminosity for the black hole 
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For M < 


, we have n max = %-tt and find 
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From eqn. (11401) for small Km p and eqn. (1141)1 for large nm p , we can conclude that the 
coefficients C impacts the black hole’s luminosity only in the intermediate range of K,m p 
noticeably. 

For the intermediate range, FIG. [l] plots the luminosity L against (■ Anm p ) 1 , which be¬ 
comes for a 2D Schwarzschild black hole with the mass M. In FIG. [H we have a = 0.9 
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and M — for w max . We plot L vs (4Km p ) _1 in FIG. H] for the usual case with C — 0 
(red line), the ones with rjC = 10 and 1000 (solid and dashed blue lines, respectively) and 
the ones with r)C = —10 and —1000 (solid and dashed brown lines, respectively). For the 
rjC < 0 cases, there are ” weird” peaks in FIG. [U which are due to the transition from 
to in Umax- However, such transitions is barely seen for the rjC > 0 cases. In our 

calculations, the luminosities are determined not only by the modified Hawking temperature 
but also the range of integration of u in eqn. (11391) . When M > - 7 ==, in the rjC < 0/rjC > 0 


. In¬ 
cases the ranges of integration are less than that in the usual case, which tends to decrease 

the luminosity. In the rjC > 0 cases, it shows from eqn. (11381) that the modified Hawking 

temperatures are lower than that in the usual case. Thus, the luminosities L become smaller 

due to the decreased temperature and the shrinked range. From eqn. (I138jl , we find that 

the modified Hawking temperatures in the r)C < 0 cases are higher than that in the usual 

case. Thus, the competition between the increased temperature and the shrinked range 

determines the luminosity. The effect of the increased temperature dominates over that of 

the shrinked range for r/C = —10 and vice versa for r)C = —1000. 

To see how the luminosities L depend on values of a, we plot L vs (4 Km p ) 1 in FIG. 
[2] for the usual case with C = 0 (red line), the ones with rjC = 10 (blue lines) and the 
ones with rjC = —10 (brown lines) with a = 0.5 (solid lines), 0.9 (dashed lines) and 0.95 
(dotted lines). Note that a parameterizes the unknown quantum gravity ultraviolet cutoff 
A. I 11 FIG. El it suggests that the r/C < 0 cases are highly sensitive to the physics at high 


energies while the r)C > 0 ones are not. I 
are super-/sub-luminal. The author in 


n > 0, rjC < 0/r/C > 0 implies that the particles 


70] has shown that the Hawking radiation with 


sub-luminal dispersion was not sensitive to Lorentz violation at high energies due to the 
’’mode conversion”. However, the outgoing black hole modes with super-luminal dispersion 
emanated from some unknown quantum gravity processes. 


VII. DISCUSSION AND CONCLUSION 

We used the Hamilton-Jacobi method to investigate the dispersive held theory models 
in this paper. Our results suggest that the thermal spectrum of radiations near horizon is 
robust. In fact, if the difference between the modified dispersion relation and the relativistic 
one was suppressed by the fundamental energy scale m p , we found that the deviation of 
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FIG. 1: The luminosity L of a 2D black hole against (4 nm p ) 1 with a = 0.9. 


the effective Hawking temperature from the standard one was also suppressed by m p . For 

TTlP 1 TTlP 1 

a particle with the typical energy u ~ -^, the deviation was given by powers of jfe. Nev¬ 
ertheless, there are some potential corrections to the effective Hawking temperature which 
are not included in our calculations: 

TTlP 1 

(a) Back-reaction effects which occurs at order fj. For a particle with uj ~ -jf, they are 
m 2 

order of jfi. However, the Hamilton-Jacobi method is incapable of computing them 
since the metric is fixed in this method. On the other hand, back-reaction appears in 


the the null geodesic method 


20 


2 l| to ensure energy conservation during the emission 


of a particle via tunneling through the horizon. These corrections lead to non-thermal 
corrections to the black-hole radiation spectrum. Note that there are some attempts 
to incorporate back-reaction effects into the the Hamilton-Jacobi method using the 


rainbow metric 


0 . 


72] 


(6) Higher order WKB corrections. In the Hamilton-Jacobi method, we take the semi- 

classical limit h —> 0 and keep only leading order terms to calculate the Hawking 

temperature. Therefore, one may wonder if the Hawking temperature could receive 

higher order corrections in h beyond the semiclassical one. The corrections has been 

2 

estimated in [73] and was given by powers of jfi. However for the non-dispersive 


34 




















FIG. 2: The luminosity L of a 2D black hole against (4 nm p ) 1 with a = 0.5, 0.9 and 0.95. 


case, several authors 74-176] argued that the tunneling method yielded no higher-order 
corrections to the Hawking temperature. Whether such arguments also work for the 
dispersive models needs to be checked. 


In this paper, we used the Hamilton-Jacobi method to calculate tunneling rates of radi¬ 
ations across the horizon and the effective Hawking temperatures in the dispersive models 
with the static preferred frame. After the spectrum of radiations near the horizon was ob¬ 
tained, the thermal entropy of radiations near the horizon and the luminosity of the black 
hole were computed. Our main results are as follows: 

• In section [TT]and the appendix, we used heuristic arguments and effective held theories, 
respectively to derive the deformed Hamilton-Jacobi equations in the dispersive models 
with the static preferred frame. Note that these methods can easily be generalized to 
any preferred frame. 

• In section QJ] the deformed Hamilton-Jacobi equations was solved for d r I and the 
imaginary part of / was obtained by computing the residue of d r I at r = 77 r . The 
assumption for our calculation was also given, which required that the singularity 
structure of d r I except the order of the pole at r = do not change after the MDR 
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was introduced. The corrections to the Hawking temperature were calculated for 
massive and charged particles to O ( m~ 2 ) and neutral and massless particles to all 
orders, respectively. It was found that corrections were suppressed by m p . 

• In section IIVI the average number and entropy for a mode were calculated for bosons 
and fermions. They could be obtained from those in the non-dispersive case by re¬ 
placing the standard Hawking temperature with the modified one. 

• In section [Vj we used the brick wall model to compute the thermal entropy of a massless 
scalar field near the horizon in UV finite and perturbative cases. In the UV finite case, 
the entropy was always finite as one approached the horizon and hence the wall near 
the horizon was not needed. In the perturbative case, a wall was put at r = ry + r E 
to regulate the UV divergence. We assumed the proper distance between the horizon 
and the wall is order of rn p . Thus, the entropies near the horizon in both cases were 
given in eqn. (1 11 7|) . We found that the leading term of the entropy depended on the 
MDR effects while the subleading logarithmic term did not. 

• In section IVII we calculated luminosities of a 4D spherically symmetric black hole 
with the mass M rn p and a 2D one. We used the geometric optics approximation 
to estimate the effects of scattering off the background. However, as discussed in 
section m such approximation might not be robust enough and hence the asymptotic 
spectrum could be significantly changed. 

Acknowledgments 

We are grateful to Houwen Wu and Zheng Sun for useful discussions. This work is sup¬ 
ported in part by NSFC (Grant No. 11005016, 11175039 and 11375121) and the Fundamental 
Research Funds for the Central Universities. 

Appendix A: Effective Field Theory and Deformed Hamilton-Jacobi Equation 

As discussed in the introduction, various approaches to the quantum-gravity problem 
could lead to the existence of MDRs. To have a MDR, one has to break or modify the 
global Lorentz symmetry in the classical limit of the quantum gravity. There are several 
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possibilities for breaking or modifying the Lorentz symmetry, one of which is that Lorentz 
invariance is spontaneously broken by extra tensor fields taking on vacuum expectation 
values. The most conservative approach for a framework in which to describe MDR is the 
effective field theory (EFT), where modifications to the dispersion relation can be described 
by the higher dimensional operators. Since we are only interested in modifications to the 
dispersion relation of the particles, we limit ourselves to the kinetic terms and neglect self¬ 
interacting effective operators when constructing the effective field theory. We also assume 
that the effective theory respects U (1) gauge invariance of the charged black hole. The EFT 
framework can easily incorporate MDR via the introduction of extra tensors. To construct 
the minimal EFT in curved spacetime, we suppose that the action of the EFT contains the 
usual minimal gravitational couplings and the EFT coefficients are constants in the local 
frame 77]. 


1. Scalar Field 


We work with a complex scalar field 0 with the mass m and the charge q. Following 
guidelines we put forth, we find the effective Lagrangian for 0 incorporating MDR can be 
written as 
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where D fl = V^ + ^Tl^, is the covariant derivative of the background spacetime, A^ is the 
electromagnetic potential, j runs over all independent operators of a given dimension, B is 
a dimensionless function of — with B (0) = 0, and, C u and Cl ,, are dimensionless extra 
tensors depending on with C0 (0) = Cj iu (0) = 0. The deformed Klein-Gordon equation 
is 
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With rotational symmetry, all extra tensors become reducible to products of a vector field 
tA, which describes the preferred frame and = 1. Thus, the extra tensors become 
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c, = c — < 


m r 


■pn ^*2fc-l^*2fc U,J " i 2k+l U ^in1 

L P, 


(A3) 


where < 7 ^ is the metric of the background spacetime, C and are dimensionless functions 
of ^-, j = (k,C) , 2 k < n, and C denotes any possible permutations of ( 1 , • • • , n), namely 
(ii,--- , i n ). To obtain the Hamilton-Jacobi equation, we make the WKB ansatz for (j) 


il 
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Defining 

f = -u» (d^I + qA „) , X 2 = T 1 - {d^I + qA^f , (A5) 

and plugging eqns. (1A4D and (1A3D into eqn. (L42lh one expands eqn. (1A2D in powers of h 
and finds to the lowest order 
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Solving eqn. 


for X 2 with respect to T gives the deformed Hamilton-Jacobi equation 
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•P/ V'f/ V"P/ n>3 

where a,/3 ,7 are dimensionless functions of with a (0) = /3 (0) = 1 and 7 (0) = 0 which 
can be determined by the coefficients B, C and C J n in eqn. (1A6IK In flat spacetime with 
A^ = 0, the dispersion relation for the scalar field can be found by inserting the positive 
energy ansatz (f> = exp ip ^ j into eqn. (L42|) . The resulting equation for p M is actually 
eqn. (1A7|) with T = tAp^ and X 2 = —p^ + T 2 , which is exact for flat spacetime with 
A m = 0. Identifying E = p = T and p 2 = —P/^p^’ M — —p^jA-t- T 2 = X 2 , we can produce 
the MDR for the scalar, eqn. ((5]) , in flat spacetime. On the other hand, the vector field tA 
is chosen to be ( A— . 0 I in curved spacetime with the metric ([6]) and the electromagnetic 

V VAH J 

potential A In this case, T and X 2 become T and X 2 in eqn. (fT3lh Thus, in the black 
hole background spacetime, the corresponding deformed Hamilton-Jacobi equation for the 
scalar field incorporating the MDR, eqn. ([5]) , is given by eqn. o. 
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2. Fermionic Field 


In the background spacetime with the metric g /JV and the electromagnetic potential A^, 
the effective Lagrangian for a spin-1/2 fermion i/j with the mass m and the charge q incor¬ 
porating the MDR can be written as 


c{f f = $ [wi'f ^) 'h 
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(AS) 
( 0 ) = 

C'/ (0) = 0, j runs over all independent operators of a given dimension, Df t = d^ + Q^ + ^A^, 
^M = 2 UJ i? b ' Bl a b, S a fe is the Lorentz spinor generator, is the spin connection and { 7 ^, 7 ^} = 
2g txv . The Greek indices are raised and lowered by the curved metric g lw , while the Latin 
indices are governed by the flat metric ?/ a &- The deformed Dirac equation is 
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ra r 


* s (7 

n>k>l,j 


2-1 — 
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k -1 


1 . . . 1 . . . _ Q 


(A9) 


With rotational symmetry, the extra tensors become 


■^Ml’”Mn 1^12 ^^*2k — l^*2k ^^2fc+l 5 


(A10) 


^Ml'"Mn ^ ^ J " i 2k-l^ i 2k^^ i 2k+l 

where £>/ and G/ are dimensionless functions of 77 , J = (fc,C), 2 fc < n, and C denotes any 
possible permutations of (1, • • • , n), namely (ii, • • • , i n ). To obtain the deformed Hamilton- 
Jacobi equation, the ansatz for -0 is assumed as 


= exp v, 

where v is a slowly varying spinor amplitude. Substituting eqn. (h411l) into eqn. 
find to the lowest order of h 


(All) 


we 


A" m 7 m v + mv + m ^ '’'^m^m^m* 

n>0,k<^,C 




Mn-l 


2fc+l 


E 


C j 


m k— l^MilMi2 9fJ,i 2 i_ iM: 


U u . • • • A" Mfc 7 Mfc+1 • • • 7 ^"V = 0, 

■2Z Mj 2 / + i Mm 1 1 


(A12) 
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where X„ = dj + qA Using {X^f = (X„) 2 = T 2 - X\ {urf) (X^) = = -T, 

and (tip 7 m ) 2 = 1, one could rewrite eqn. (\A12\i as 


m 


m r 


Xu + rnC] — ) Uu + h 2 (T, X 2 ) Uu + g (T, X 2 ) X 


7 m v 


m + mC'2 



+ />, (f,A' 2 ) 

(A13) 


where C z are dimensionless 


functions of — with C\ (0) = 


E 

2p+q>l 


2p+q-l 


and g t T, X 2 = 


E 

2p+q>0 


jP,q ( at ^Tt^T 2 X 2 ) 
,^P+9 




0 , hi (T , X 2 ^j = 

The coefficients C t , 


h™ and g p,q are determined by B 3 n and C J n from eqn. (IA12D . However, the detailed relations 
between them are irrelevant here. Multiplying both sides of eqn. (1A13j) from the left by 
(X )t + mCiU^ + h 2 u fl + gX^) 7 ^ and then using eqn. (IA13D and { 7 ^, 7 v } = 2 g^ v to simplify 
the RHS, one gets 


v, 


(X p + mCiu p + h 2 u fl + gX p ) (X^ + mCiu^ + h 2 u ll + gX^) v = (m + mC 2 + hi) 2 v. 

(A14) 

Since v is nonzero, eqn. (L414I) gives 



where 


Z TD \ m 

■m B I — 
rn„ 


■7nc(—)f+ V 

V P/ n>2,fc<f,C 


C 2 , 

n \ m z 


T 2 — X 2 ) T 


pn—2k 


771, 


n—2 


- = 0 , 

(A15) 


B = -Cl + C 2 2 + 2C 2 , 

C = 2C\ + 2g°'°C 1 + 2Cih°’ 1 - 2 (1 + C 2 ) h 0 {\ 

CJj f hp2 r pn—2k 

- ^X 2 -= 9 (t 2 - X 2 ^j + 2 (h 2 + gh 2 ) f - h 2 + h 2 

+2 771 ( g - g°’°) Cif + 2ttiCi (h 2 - h^T) - 2m (1 + C 2 ) (h t - . 

It is noted that the form of eqn. (1A14D is the same as that of eqn. (L46j) . Thus, the argument 
and result below eqn. (jAL7j) can also apply to a spin-1/2 fermion held. 
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